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The formalism for Schwinger-Dyson equations of Wilson loops, as developed by Eguchi, 
Weingarten, and Foerster, is extended to include connected correlators. These should allow 
one to study propagators and interactions of glueball fields without having to handle gauge 
dependent fields such as gluons an ghosts. Some steps towards a solution are made. 



p ,! 1 Introduction 

CD ■ 

i-Q I Schwinger-Dyson equations are a powerful analytic tool to study quantum field theories beyond 

perturbation theory. In writing down the equations, no approximations need to be made, which 

^SJ I means that the equations contain, in principle, all physics — both perturbative and nonperturba- 

<*" • tive — exactly. 

^^ I In a general quantum field theory, the equations follow from Stokes' theorem applied to the 

path integral: 

q: /[d^]^e-^[^i+^-^ = o. (1) 

^^ ■ From this master equation, the entire tower of Schwinger-Dyson equations can be written down 

Cn I by expansion in the source field J{x) and equating the expansion coefficients with zero. This 

leads to equations involving full (i.e. containing disconnected parts) n-point functions, and with 
some further manipulations they can be reduced to equations involving only connected or only 
IPI Green's functions. 

In the case of Yang-Mills theory (and, by extension, of QCD) the role of the field (p{x) is 
CTJ ■ played by the gluon field AJJ(a;). This means that the Schwinger-Dyson equations obtained in this 

way will contain gluon propagators and gluon interaction vertices. But what, in the end, interests 
us most are not so much gluons but rather bound-states like glueballs (and, in QCD, hadrons), 
meaning that after solving the full tower of equations — supposing that we are smart enough to 
do so — we would still be facing the formidable task of putting all elements together to extract 
the physical quantities of interest. One can therefore wonder whether it would not be possible to 
write down a set of equations involving only glueballs, thus sidestepping at least one leg in the 
quest to a better understanding of Yang-Mills theory. 

This leads us to the question of the gauge symmetry. The gluon field has a nontrivial gauge 
transform, which makes it rather tough to get rid of the gauge dependence when starting from 
something like ^. Ideally, one would like to take the functional derivative in some gauge invariant 
way, such that only gauge invariant quantities show up in the resulting equations. These gauge 
invariant operators would then be interpreted as glueball operators, and the goal mentioned above 
would be reached. 

It seems, however, that such an ambitious program does not come easily when working in the 
continuum. Work in that direction has been done in times past [U El El H], but it seems like 
results are come by most easily when introducing a lattice regulator. This means that manifest 
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Lorentz invariance will have to be sacrificed, to be subsequently recovered when we are able to 
take the continuum limit. All work done in that direction in the past (see, for example, [51 [SI [71 |H] ) 
focused on computing expectation values of Wilson loops, hoping to get a better understanding 
of confinement. For a review see 9^ and references therein. Section [2] reviews the formalism of 
Schwinger-Dyson equations for Wilson loops on the lattice. 

In my work I aim for connected correlators, which should give more insight in the spectrum of 
the theory and in scattering amplitudes. In section[3]the equations for different connected and IPI 
n-point functions are written down and some issues are discussed. In section [4] I discuss different 
approaches aimed at solving the equations at three-level, i.e. in the large- iV limit. Finally, section 
[5] concludes this paper with an outlook at future possible avenues of research. 

2 Schwinger— Dyson on the lattice 

2.1 The lattice functional derivative 

In order to write down Schwinger-Dyson equations, one needs to define a derivation operator 
Stokes' theorem can be applied to. On the lattice, the integrations are integrations over the gauge 
group, and a derivative exists having the desired properties J 101: 
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where r° are the generators of the gauge group under question, obeying the commutation relations 
[t'^,t^] = 2if'^^'^T'^ with /"'"^ the structure constants, and the dots denote other stuff the quantity 
Q may depend on. Mark that the derivative carries a gauge group index a. When computing 
the integral of such a derivative, the result will be zero due to the invariance of the integration 
measure under U — > VU with V an arbitrary gauge group element. A right derivative can be 
defined in a completely analogous way. 

In practice, this functional derivative operates as 
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and it furthermore obeys Leibniz's rulelj When working on a Wilson loop containing the link 
matrix Ux,fj,, the derivative will "cut open" the loop and introduce a group generator. One can 
depict this as: 
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Here, the arrow shows the orientation of the Wilson loop (i.e.: the loop in the first identity contains 
the matrix Ux,fi, the one in the second identity contains U^ ^). The dot with index a has been 
added at the place where the group generator iT"/2 has been inserted. A derivative 5a;+/i,-^p 
working on the same link matrix but coming from the other side, gives an analogous result but 
with the signs flipped around. 

In order to find a gauge invariant result from the lattice functional derivative, one needs a 
way to close the Wilson loops that have been cut open by the derivative. Given the well-known 
identity for the generators of the gauge group 



E^"^^" 



Ta/^T^S = '^SasSfj^ - —SafsS^S 



(5) 



^The derivative operator itself commutes with all other quantities. 

^I will use the notation fj, for the vector with length the lattice spacing in the direction fi throughout. 



where N is the number of colors (in S\J{N)), it is clear that having two generator insertions at 
the same lattice site and summing over the color index will close the loops: 
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One sees that, in order to have closed loops as a result, it suffices to already have some generator 
with same gauge index a present at the lattice site where the derivative is going to be taken, after 
which one has but to sum over the color index. 

When the derivative in this gauge invariant combination hits the piece of Wilson loop where 
the generator matrix was inserted, another useful identity can be used: 



E 



»t" N'^ - 1 
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This is an immediate result of the identity ([5]) . 



2.2 Towards the equations 

All this together bring us to the strategy to follow for writing down Schwinger-Dyson equations. 
The master equation iso 



Y. h^UKAQ^^U]n[U]e- 
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Here, Qx[f^] is some Wilson loop that has a gauge group generator insertion at the site x, TZ\U] is 
a product of some Wilson loops (or none), and S\U] is the Yang-Mills action. In this paper, the 
action will be taken to be the ordinary Wilson action: 



S = 



N ^ 
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TV 



PA^) 



(9) 



where P^^{x) is the one-plaquette Wilson loop in the (/t, v) plane centered on x. 

A first term encountered when working out ([5]) is the one where the derivative hits Qa;[t^] 
in the link where the gauge group generator insertion has been placed. The identity ([7]) will be 
triggered, and the resulting term is nothing but 



N' 
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A second contribution comes from the derivative hitting the loop Q in another instance of the Ux,fi 
matrix, as can happen if the loop doubles back on itself. A first term from this is the one coming 
from the first term in ([6]), where the generators are just deleted and a contribution proportional to 
(|10p arises. An additional term comes from the switching around of the two pieces of loop, which 
will result in the loop Q falling apart into two sub-loops. As an example, consider a Wilson loop 
winding twice around an elementary plaquette, with one insertion of a generator and a derivative 
at the same site: 
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^In the typical way of constructing Schwinger-Dyson equations, one introduces a source term, and different 
equations are obtained by taking derivative with respect to the source. In the case at hand it seems more convenient 
to put the relevant fields already in the path integral instead, as one of the loops should contain a gauge group 
generator. 
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Figure 1: The one-plaquette operator P and the five two-plaquette operators: the double plaquette D, 
the seat 5", the rectangle R, the flipped rectangle F and the folding chair C. 

Here, the a;-axis is taken to point to the right. 

If, in the master equation ([5]), one of the loops in TZ[U] contains the link the derivative is 
working on, there will be additional terms from this. A last contribution comes from elementary 
plaquette Wilson loops coming from the action in the exponential. These will give higher n-point 
functions containing these elementary plaquettes as well as terms wherein the loop Q has been 
augmented with these plaquettes. 



3 Constructing the equations 

3.1 One-point functions 

As a first case, let us consider the equations for vacuum expectation values of Wilson loops. These 
have already been derived in the literature, but for completeness I will repeat the derivation using 
the notations introduced in this paper. First let me introduce a slightly deviant convention which 
will simplify the look of the equations later on. With P^i, (x) I will denote the Wilson plaquette 
centered at x turning in positive direction in the coordinate system (/t, 0) . This is depicted in 
Figure Hal The vacuum expectation value of the plaquette divided by the number of colors will 
be denoted P. Mark that an operator like P,^fj.{x) is equal to P*^,{x). 

There are some more operators needed in this paragraph. They are also drawn in Figure 
[TJ For these operators I will again use the convention of indexless letters to denote the vacuum 
expectation values, thus: D = {jj-Df^„{x)) etc. Furthermore, it will now and then be useful to 
allow indices to take negative values, with which I mean that the direction in question is taken in 
opposite sense. This means that, for example, P-^i,{x) — P*{x). 

3.1.1 The elementary plaquette 

As a first example, let us consider the master formula ([8]) with Q[U] the elementary plaquette 
Wilson loop P centered at the origin, TV[U] absent, and the derivative taken along one of the links 
of the loop under consideration. For concreteness, let the loop be in the (x, y) plane, with the 
derivative working on the site at —^{x-\-y) along the x-axis. The starting point is: 



Y^j[du]5i{-^:^)(p-^ 
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The superscript at the plaquette indicates that it has been cut open with a group generator at 
that lattice site. 

The first term in the equation is the one hke in (flUl) : 

7V2-1 

where one factor of N has been absorbed in the vacuum expectation value of the loop. 

The second contribution comes from elementary loops coming down from the action. Firstly, 
there is the term where the loops are not switched around (i.e. the contributions from the first 
term in ([HI)), which introduces all elementary loops containing the link of the functional derivative, 
with a minus sign if the loop flows in opposite direction from the derivative: 

^(^..(0) E {P^.i^) - P,A^))) ■ (14) 

^ \ ±fi^x I 

The sum over /i is here over both positive and negative senses, in all directions except along the 
X axis. If now we introduce the propagator notation 

A^.,«A(a: -v)^ (P^.(x)Pak(2/))co„„ - {Pt.Ax)Px.{y)) - N^P^ , (15) 

the above result can be rewritten as 

-1^ E (a..,..(^)-A..,..(^)) • (16) 

Then there is the part where the loops coming from the action and the original loop put in the 
path integral are linked with each other, through the second term in ([5]), giving 



N 
2^ 



^^iD-l + F-R + 2{d~2)iC-S)) (17) 



where d is the spacetime dimensionality. The term of N comes from a tr 1, caused by a loop being 
reconnected to its mirror image. 

All together we find the following equation: 

N'^ -1 1 

7^P + TY^i^D -Ai+Af-Ar + 2{d - 2)(Ac - As)) 

2 Zg''I\ 



N 
2^ 



^iD-l + F-R + 2{d-2)iC-S)) = 0, (18) 



where Aq is a shorthand for the propagator between two plaquettes that share at least one link, 
and which are positioned such that linking them with each other would give the operator O, as 
for example A^i = Axy.xy{0). This gives a relation between the values of certain condensates and 
the propagator of the operator P^^ at infinitesimal (lattice-spacing) distance. 

In order to make contact with the large- A^ expansion, let us introduce the 't Hooft coupling 
A ~ g^N . Using this, the above equation can be rewritten as 

^- ^) P +\iD + F - R + 2{d-2){C - S)) ^^ 



1 1 

772 A 



2 JAz5 - Ai + Af - Afl + 2{d - 2){Ac - As)) . (19) 



This equation can be schematically depicted as 

• + ^=#= • +]^0 • (20) 

The left-hand side contains the one-plaquette condensate (single line) together with several other 
condensates (double line), which are the unknowns of the equation. At the right-hand side one 
has a tree-level value for the condensates, and a loop-correction which is suppressed by a coupling 
constant equal to 1/N'^. Mark that the first coefficient at the left-hand side receives a 1/iV^ 
correction which does not originate from any loop-like diagram, but which is rather similar to a 
finite rcnormalization. 

3.1.2 Other Wilson loops 

Equations for other one-point functions can be written down in a completely analogous way. 
Several issues immediately arise. 

When writing down the equations for, for example, the double plaquette D condensate (see 
Figure [Tb]). it will contain a nonlinear term due to the recombination coming from the second 
contribution shown in ()11|) . The presence of such nonlinearities in the condensate equations 
seriously complicate solving them. 

Another point of note is the fact that there exists a certain arbitrariness in the construction of 
the equations. In order to write down the equation for a certain condensate, a site and link in the 
relevant Wilson loop have to be chosen, and a different choice will, of course, lead to a different 
equation. The solution must still be the same, which it can be verified to be in those cases where 
explicit expressions for (parts of) the solution can be obtained, like in the limit of large 't Hooft 
coupling A within the large- A'^ limit. In the following I will systematically symmetrize the choice 
of site and link over the entire loop under consideration, so as to get more symmetric equations. 
What extra information, if any, can be gained from making different choices is still to be explored. 

3.2 Two-point functions 

Let us now consider the Schwinger-Dyson equations obtained by putting two one-plaquette op- 
erators into the path integral. One of these plaquettes is positioned at the origin as in the one- 
plaquette example, and the second one has an arbitrary position x. The expression we consider 
is 

i E E l\-'^U]5t {p:j{Q)P,s{x)e~') = . (21) 

links ; a •' 

The j'^i is a symmetrization over the four links of the square plaquette at the origin, which will 
make the equations manifestly symmetric under rotation over a right angleo 

The equations can be straightforwardly written down following the same procedure as in the 
previous section. The main difference is the extra contribution arising whenever the second pla- 
quette Prys{x) has some link in common with the first one. In that case, two extra term appear due 
to a functional derivative working on this second plaquette, giving one term where the plaquettes 
are left unchanged and one term with the two Wilson loops merged into one. Schematically, the 
result looks likeH 



^(-)("+f+i.0)4-(^_^-(22) 



All n-point functions with blobs denote full n-point function, which contain the disconnected 
parts. In the above equation, 6{x) depicts schematically the Kronecker deltas coming from all 



^Unsymmetrized equations must also be fulfilled by an eventual solution of the theory, but, pending the de- 
velopment of the necessary techniques to find such a solution, it seems one gets most physical intuition from the 
symmetrized equations. 

^The full expression is quite lenghty, especially for general indices and positions, but the interested reader should 
not encounter too many difficulties in deriving it for himself. 



different configurations where the two plaquettes touch each other. The constant a is one if the 
two plaquettes in the equation exactly cover each other but have opposite orientation, and zero 
otherwise. 

The next step is to separate all n-point in a connected and a disconnected part. The discon- 
nected parts coming from all terms — except the one with the Kronecker delta — cancel due to 
the equation for the one-point function. Next, we go to Fourier space using the transforms 

Kv) - E^"''/(^) ' /(^) = /7''"''^~^^^(S^ ' ^^^^ 

where the integral in the second identity is over a hypercubic box of size 27r. The result is a 
system of linear equations for the propagators Aa^_-y5(p). It is possible to decouple these equations 
according to charge conjugation by taking symmetric and antisymmetric linear combinations of 
the Lorentz indices. Defining P/^ — \{Ppiv ± -Pi^a*) (which are nothing but the real and imaginary 
parts of the plaquette operator P;^,y), one gets the new propagators 

A^l,„^(x) ^ (P^^(0)P*^(x))_ . (24) 

With these definitions, there are four propagators, but the ones mixing different charge conjugation 
states must be zero due to the charge conjugation symmetry. This turns out to be consistent with 
the equations. For the two nonzero propagators we get: 
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where the following shorthands have been used: 

, I , , (Jq-v cos ^- cos ^ + (5fl-, cos ^ cos ^ + ^^^ cos ^ cos ^ + (5«5 cos ^ cos ^ 
C/3;75(P) = -^ —^ (26a) 

(5q7 sin ^ sin ^ - (5^^ sin ^ sin ^ - (5q5 sin ^ sin ^ + (5/35 sin ^ sin ^ 

and where furthermore the operator Q^ui^x) has been defined as four times the Df^^{x) operator, 
plus the sum of all F and C operators containing Ppj,(a;), minus all the R and S operator containing 
this same plaquette. Furthermore, the connected three-point function has been introduced as 

Tap,js,edx,y,z) = {Pap{x)Pys{y)P,^{z))^^^^^ , (27) 

and with similar expressions for symmetrized and antisymmetrized three-point functions. Its 
Fourier transform has been defined as 



x,V 



The above equations can be depicted as 



o- 



where now the propagators are connected and the three-point vertex is the IPI one. The matrices X 
and Y have been introduced and depict the expressions that can be found in the explicit equations 
above. These equations can be compared to the standard form of Schwinger-Dyson equations. 

The matrix X would ordinarily be the bare Euclidean propagator m^ + p'^ , eventually com- 
plicated with some Lorentz index structure. This is quite similar to what we have here for the 
charge conjugation odd sector in equation (|25bl) . with the main difference that, on the lattice, 
momenta are replaced by sines of momenta. The charge conjugation even sector has, at this point, 
no momentum dependence in its bare propagator yet, but the inclusion of more complex Wilson 
loops beyond the elementary plaquette immediately remedies thatO The second term in the above 
equation depicts mixing with other fields. In order to extract the physics in the large- TV limit one 
should diagonalize the full propagator matrix, containing all mixings of all different shapes of Wil- 
son loops. The Y tensor would, in ordinary Schwinger-Dyson equations, be a Dirac delta function 
in position space, or a constant when going to Fourier space. In our case, it has some nontrivial 
momentum dependence due to the mismatch between which fields are in the path integral (the 
link matrices) and which fields we are considering n-point functions of (Wilson loops built from 
link matrices) . The last term in the equations is a loop correction due to a three-point interaction, 
governed by a coupling strength 1/A^. Naively one might, at this point, be led to believe that the 
effective theory describing glueballs will only have a three-point vertex. Explicit consideration of 
higher n-point functions (for which see the following section) shows, however, that matters are 
more complicated. This can already be forseen due to the l/N"^ corrections the Y matrix receives 
(see the equations (I25p ). which superficially look like a tadpole correction from a four-point vertex, 
but which, all things considered, do not have the correct form to be exactly that. 

3.3 Higher n-point functions 

3.3.1 Notational conventions 

In order to discuss higher n-point functions in a more structured way, I will introduce some matri- 
cial notation. Let us consider the Schwinger-Dyson equation for the n-point function containing 



®Also, the momentum dependent term in the charge conjugation odd sector has the wrong sign, which also 
drastically changes when including the mixing with other Wilson loop shapes. 



the Wilson loops A{0), B{x), C{y), etc. The functional derivative is taken at links of the loop A, 
which is cut upen by a gauge generator as necessary. 

First let us suppose we fixed some choice determining how the functional derivatives are going 
to be taken. This consists of choosing the links of A we will use, and a linear combination to 
sum the different contributions. Let us furthermore suppose that the coefficients in the linear 
combination add up to —2/N at leading order in 1/A^|3 so chosen to have less clutter later on. 
Let us now consider the different terms appearing in the path integral. 

First there are the different ways the derivative can work on A{0). This will lead to contribu- 
tions containing just ^(0), and contributions where the loop A falls apart into two different loops, 
analogous to what is shown in (|lip . Whichever choice one made concerning the links the functional 
derivative works on, the coefficient of the contribution containing just A(0) will always be one at 
leading order for large N, due to the scale of the linear combination chosen|f| To describe the 
contributions from the loop A falling apart into two loops J and K, let us introduce the notation 
a as: 

^ Y^ aA(o),j{3),K{k)J{i)K{k) . (30) 

J(3).K{k) 

The symbol a may, of course, depend nontrivially on iV, but with the definition above its leading 
behavior at large N will be of order N^ . The symbol a can always be defined to be symmetric in 
its two last indices, but this is not strictly necessary. 

If one of the other Wilson loops under consideration S, C, etc. is well-positioned, the functional 
derivatives may work on it as well and give two more kinds of terms. The first kind is where the 
two loops A and, for definiteness say, B are left as they are. These will be written af|j 

- ^A(0)i?(x)/3^(o)'^(-) , (31) 

where /3 is again order iV° at large N , and where the sign is chosen to have /3 be positive if the 
links of A and B hit by the functional derivative go in the same direction. The second kind are 
those where the loops A and B are linked with each other to form one larger loop J: 

j^YlA(a),B(x).,j(j)J{i) ■ (32) 

It is possible for A and B to cancel each other and give the unit operator, as happens whenever 
they follow the same loop at the same spacetime point but have opposite orientation. This will 
not be included in the 7 symbol and will instead be treated separately. 

Finally there will be contributions from the functional derivative working on the action, bring- 
ing down elementary plaquettes. These can be left unlinked with A: 

- ]^^(0) E 4S^b) > (33) 

p(p) 

or linked with A: 

^ECa(o)„70)J(j)- (34) 

JU) 

The notation ^A(a;).s(y) will b^ used as a Kronecker delta for Wilson loops, and it is one if A 
and B are identical loops taken in the same direction, and at the same spacetime point — and 



^If this is not the case, one has but to multiply the coefficients with a suitable number to adjust it. This will 
only affect the equations in multiplying them with a constant. Only when a certain choice leads to a zero sum of 
coefficients will one cncoutcr a problem in doing so. I will in the present work ignore this degenerate case. 

*Subleading terms of this coefficient can always be made to disappear by adjusting the linear combination. 
The difference between putting the indices above or below has no specific meaning, but is chosen such that a 
semblance of the Einstein summation convention can be applied to indices put below, if one interprets condensates 
to have a lower index. 



zero otherwise. Furthermore propagators and vertices will from now on be written in matrix-like 
notation, as l::^A(x),B(y) and TA(x),B(y),c(z)- 

As an example, the equation for the one-point function of the Wilson loop A will now be 
written as 



JU),K(k) JU) 



+ ]V2 ^ "A(o),JO%K(fc)Aj(j),K(fc) - j^ Yl ^P{p)^M0),P{p) - ■ (35) 

J(j),K(k) Pip) 

For the propagator equation, and also for further use, it comes in handy to define the following 
matrices: 

'^A{x),B{y) ^SA(x),B{y) +YaA{x),.J{j),B(y)J+ ^ aA(x),B{y),K{k)K 
JU) K(k) 

- J^Yl ^P(p)^Pip)My) + -y^U{x),B{v) ' (36a) 

P(p) 

'^A(x),B{y) = - ABP'^'^'^^'-^'^y'' + ^ 'yA{x),B(y),J(j)J ~ SA{x),B'{y) , (36b) 

JU) 

and the propagator equation becomes very simply: 
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2J '>^A(Q)J(j)^J(j)^Bix) + yA(0),B(x) + -^ 22 ^A{0).,J{j),K{k)'^''j{j)^K{k)M(x) 
J if) J(j).K{k) 

~ iVA ^ ep(p)r^°o) P(p) B(^) - ^A^(o),B(x)/5'^^°^''^^^^ = . (37) 

Pip) 

This means that the inverse propagator at tree level (i.e. in the large- A^ limit) is equal to — Y^^X. 
Mark that Y has to be invertible for this to exist. It can easily be seen to be so in the limit of large 
A, as the matrix reduces to —SA{x).B*(y)- I will in the following always suppose it to be invertible 
for low A as well, although I have no formal proof of this. 

3.3.2 Three-point functions 

The procedure to write down the equations for the connected three-point vertices is completely 
similar to the one for propagators. When one goes from connected to IPI, however, two things 
are out of the ordinary. 

As is normal in Schwinger-Dyson equations, one of the terms in the equations for the connected 
three-point functions is 

2_,^'4(0),JO)r™(""B(^)^C'(y) • (38) 

J(j) 

When writing the connected vertex in terms of the IPI vertex and propagators, one gets for this 
term: 

X-#-tr ' (39) 



-< 



Then, the product XA, together with all loop corrections that are not IPI, can be replaced by 
their value from the propagator equation. In the ordinary Schwinger-Dyson formalism, that would 
be a delta function, leaving just the IPI vertex with two propagators sticking out. In the case at 
hand, a Y matrix appears instead. 

10 



A second divergence from the ordinary Schwinger-Dyson formalism consists in deleting the 
propagators sticking out to the right in (1391) . Most terms in the three-point equation have exactly 
the same two propagators sticking out, as is ordinarily the case, such that those just get deleted 
to come to the IPI equations. However, there are also some terms with the form 



N 



■(IX. 




(40) 



plus the same with the two rightmost legs interchanged. The circles with /3 and 7 inside depict 
the symbols defined above. These terms do not have the two propagators sticking out at the right, 
meaning that they will eventually get inverse propagators multiplying them. If we are interested 
in the large- TV values, this inverse propagator can be replaced by its tree-level value of — Y~^X. 
The final result in the large- A^ limit is 



^0) 



N 



'A(0),C(y) 



^J^A(0),JO)^-l 



J{j).B{x) 



^ Yl ^A(a),J(j),C(y)^.jlj)^B(x) 



J(J) 



N 



■'(J) 



-^SAio),B(^) E 4^SMp)My) + (B{x) o C{y)) + 0{^) . (41) 



iVA 



p(p) 



This also confirms that the three-point vertex has order 1/A^. 

3.3.3 Four-point functions and higher 

The issues encountered when writing down the equations for higher n-point functions are similar 
to the ones of the three-point vertex. When going from the equations for the connected vertices 
to the ones for the IPI vertices, one has to use the equations for the three-point vertex to get rid 
of certain terms that are not IPI, but the terms like the ones in (|40|) do not match and introduce 
some extra terms. These look like 



^% 




N -(7 




•r 



(42) 



and their permutations. Due to these terms, the tree-level value of the four-point vertex will 
depend on the three-point vertex. This will carry over to higher n-point vertices. The final result 
for the four-point vertex for large N is 



51 '^A(Q),JU)^]F{]),B{x),C(y),D{z) " - TH^ A{0) M^) P^^^^ '^ '^^^ ^ b\x) ,C {y) 



JU) 



N 



B{z),C(y)M^) ~^ J^'' 



^IPI 

JU).C{y)M{z) 



JU) 
^ ^^7A(o),J(i),B(x)i 7(j"),c(y),r>(2) + (permutations of B{x), C{y) and D{z)) + O(^) . (43) 



N 2^' 



JU) 



iipi 



The four-point vertex has order 1/N'^. Higher n-point vertices will have order l/N"^' 
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4 Towards concrete results 

In this section I present several concrete results one can get out of the equations formulated 
earlier. A full solutions is, even in the large- iV limit, not straightforward to obtain. For example, 
in order to obtain the spectrum, one should compute the matrices X and Y in p6p . after which the 
product — Y^^X (or its inverse) should be diagonalized. As both matrices are infinitely big, this 
is a rather formidable task to do with only brute force. And without having at least the inverse 
of the propagator, interaction vertices cannot even be considered. 

One regime that is readily accessible with only brute force is the limit of large 't Hooft coupling 
A. The physical limit is, however, the one of A going to zero, as this will be the continuum limitP"! 
In the large-iV limit, there appears to be a phase transition when going from large to small A, 
making a strong coupling expansion less usefull as a way to approach the physical limit [TT] . 
Furthermore, one has to go to rather high order to see any difference between different numbers of 
spacetime dimensions, which means that the rather trivial results of two-dimensional pure Yang- 
Mills theory will be found until one can gather the patience to go through the ever more complex 
higher-order computations. 

Another case that is partially accessible is the one of two spacetime dimensions. Again it is 
not physically the most useful case, but it can serve as a playground to test the formalism. 

A possible way to simplify the task of extracting physics is by doing a change of variables in 
field space. It seems to be possible to get rid of all inverses of Y in the tree-level expressions for the 
inverse propagator and for the three-point vertex, and furthermore all the matrix products that 
are left become finite sums. This means that arbitrary terms in the Lagrangian density describing 
the effective theory can be written down in finite time without the need for further artifice. This 
does not yet yield the physics we are most interested in, as the propagator will still be an infinitely 
big matrix function of the unknown one-point functions, but it is a promising avenue for further 
exploration. This is discussed in paragraph 14. 2 1 

4.1 Pure Yang-Mills in 2d 

Let us first consider the case of two spacetime dimensions. Following the usual steps (see, for 
example, |11) ) it is possible to reduce the theory to contain only one link matrix: 

Z= /[dt/]ex(f^f^+ft^-^) . (44) 



The observables of the theory are described by traced powers of this one link matrix: 

trC7". (45) 

The equations for the vacuum exp ectation values hereof in the large- A'^ limit have already been 
considered and solved in [TU [5] O 



'(l-^A)(^^ + %iii^) (A<2) 



W„ = (^trC/") = r 2 ^V "("+!) "("-!) ) ' ' , (46) 

ij<5„,i (A>2) 

where the Pn{x) are the Legendre polynomials. At the point A = 2 a third-order phase transition 
happens [TT]. The Wn are once continuously differentiable at that point. 

For A > 2, the system of equations for the propagators simplifies considerably due to most 
one-point functions being zero. After splitting the system according to charge conjugation, the 
solution can be readily found to be diagonal in both sectors: 

TTi TTli 

^mn = '77°rn,n , ^m,n — ~~;y(>mn ■ (47) 



^"As I have not spoken a word of renormalization, A is always the bare couphng, which is hnked to the lattice 
spacing through the renormalization group equations. 

^'^In this expression, the case of ra = 1 in the A < 2 region can be obtained by taking the limit n — )■ 1. The 
explicit result is Wi = 1 — ^A. 
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Interaction vertices all turn out to be zero. Mark, however, that the absence of interactions at 
tree level does not imply the absence of \/N'^ corrections, as the bare propagator, and indeed even 
the vertices, receive "finite renormalizations" , thus leading to non-trivial behavior at nonleading 
order. Indeed, explicitely going through the equations to compute the corrections uncovers, for 
example, a nonzero correction to the one-point functions. In order to compute this correction, one 
can use a strong-coupling series as an ansatz for the correction: 



11°° 1 

;^TTrE°^»^'T^' (48) 



after which a recursion formula is obtained for the coefficients a„,i. This recursion formula can be 
recognized to be the one for the Catalan triangle in slightly shifted form. The one-point functions 
are then: 

Wn ^ \5^, + ^ (^5„, + l^^l - ^) ) + O(^) . (49) 

The case of A < 2 is more complex, and it is not possible to solve the equations merely by 
squinting at them. Instead one can use the technique of the generating function [T^l [5] to rewrite 
the equations in term of the unknown generating function of the propagators ^^^^ ^ x"^y^ /S.mn- In- 
order to solve this algebraic equation, one needs the input Ai„. In the charge-conjugation even 
sector one happens to have the identity 



^tn^--—Wn, (50) 



YdA 
which follows from deriving the path integral with respect to A. The result can be found to be 



A++ = 



I >L >L — >L 22^ _ 3Af^ -L ISA"* \ 

' 8 24 82^ 32 '"'' 

a!. _ Ai 9\2_o\3 I 9Af 9A^ 33A-^ , 21A'^ 9A'' 

2 4 ZA ^A -|- ^g 2 4 "•" 4 8 

9Af, _ 3Af, I 15A"^ 9A^ _ 33A^ , 21A'' _ 9Af, 81A^ _ r,y x3 i 459A"^ _ 27A-'' , TSA" 

8 2~'~32 2 4~'~4 8 8 '•^"^16 2'" 32'' 



V ; : ■■■} 

(51) 
This connects continuously (but with discontinuous derivative) to the result for A > 2. I have not 
yet found any way to compute the charge-conjugation odd propagator. 

The charge-conjugation even connected vertices can be computed in the same way. The three- 
vertex turns out to factorize: 

rtl„+ = ^7a™7n , 7« = -^(^;(l-A) + P^i(l-A)). (52) 

It is discontinuous at A = 2. Extracting the IPI vertex would involve inverting the propagator, 
or alternatively (if one starts from the equation (J4T|) for the IPI vertex) inverting the Y matrix, 
neither which I have been able to do so far. 

This rather trivial example of two-dimensional pure Yang-Mills theory shows the soundness of 
the formalism. From the equations, the symmetry of the n-point functions under interchange of 
their legs is not at all trivial, as one leg is alway specially singled out. The results cleanly recover 
the necessary symmetry. The example also shows the difficulty in extracting physical information, 
as infinitely big matrices have to be inverted. In the higher-dimensional case, this will only be 
worse, as it will not be possible anymore to reduce the theory to as simple a set of fields as in (|45|) . 



4.2 Beginnings of an effective theory 

It turns out a certain simplification happens when performing the following change of variables in 
the effective action: 

A(x) = Y. '^B'iv),A[-)B'{y) , (53) 

B'{y) 
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where the fields with a prime are the new ones. For the IPI n-point functions this results in the 
change 

.K3),K(k) 

^'a(x).B(v),C(z) = X! ^A(x),.J(3)^B(y),K(k)^C(z),L(l)^J{3).K{k)Ml) ' (^^b) 

J(j),K(k)Ml) 

and so on for higher vertices. The new inverse propagator then becomes at leading order in 1/A^: 

^'a{o),b{x) = ^ X! ^^(o),J(j)Yb(x),jo) + 0{^) . (55) 

J[i) 

Not only is it no longer necessary to invert an infinitely big matrix, but furthermore the sum over 
Wilson loops J is finite, as X^(o),J(j) is, for any A, only nonzero for a limited number of J(j). 
With the above formula, one can, in principle, write down any propagator term from the efi^ective 
Lagrangian describing glueball dynamics. 

One could argue that the change in variables ([55)1 is, in general, highly nonlocal, thus possibly 
introducing spurious zeros into the inverse propagator, or a nonphysical particle into the theory. 
This indeed happens, as can already be verified by considering the two-dimensional theory without 
reducing it to one link (i.e. two-dimensional pure Yang-Mills theory in its most naive form). In 
the charge-conjugation even sector, the inverse propagator for the operator P{2{x) turns out to be 

{A'-^)++p,{p) = 2(1 + 3Wl + W2) - 2Wf Y. sin^ 2^ . (56) 

A' 
If this is compared with the inverse propagator of a free, massive scalar particle on the lattice 

m^ + ^sin^^, (57) 

it looks like a physical particle must be present in the theorvri while we know that the two- 
dimensional theory with only glue has an empty spectrum. In order to be able to discern between 
real particles and fictitious particles introduced by the change of variables (|53l) , one has to take 
into account the interactions between the particles, as such a fictitious particle will have no on-shell 
interactions. 

With the above change in variables, it is also possible to compute elements of the three-point 
vertex at tree level: 

1 ^^ 



^'A{Q),B{x),C{y) - ^ ZZ "A(0).JO%X(fc)Y 

./(i),-ff(fc) 



^ 2^ (^A{0)JU)M{k)'<B{x),J{j)''C{y),K(k) 
.h3),K(k) 



f^"^_^^ ' '^y^h-^uj ■ '-^yyjj^uj ■ j^t / ,- '- ■ ■'^\x),.jyj) ■ C{y),A{f)) 

^ N E ^A(0),J{j)^K{k)'>^B(x)J(j)"^Ciy),Kik) " ]^ E ^P{p)^ B{x)MO)^C{y),Pip) 
J{j),K(k) Pip) 

+ (Bix) ^ C{y)) + 0{^) . (58) 

Again, whatever be A, all sums will be finite. This allows one to write down any cubic term in 
the efi^ective Lagrangian with a finite number of computations. If one applies this to the case of 
the ^^12(2;) operator in two dimensions, one finds that the vertex is also momentum dependent, 
thus leaving open the possiblity of the vertex to have a zero whenever the inverse propagator has 
a zero. In order to verify this, however, the full inverse propagator should be known, thus still 
calling for more powerful techniques to tackle the problem. 

The formula for higher n-point functions does not undergo the same simplification. 

^■^Including mixing terms with other fields may lead to a more complex picture, but it will not make the pole in 
the propagator disappear. 

14 



5 Conclusions and outlook 

In this paper, Schwinger-Dyson equations for correlators of Wilson loops were developed. It turns 
out to be possible to write down equations for glueball propagators and interactions, and the 
coupling strength governing loop corrections is the inverse number of colors 1/7V. This means 
that, in the large- A'^ limit, the equations simplify to equations for tree-level n-point functions, and 
1/7V corrections can be systematically written down. 

In order to extract physics, however, more work will still be needed. Even at tree level the 
tower of equations is infinitely large, and only in the trivial case of two-dimensional pure Yang- 
Mills can the solution be partially extracted without too many additional difficulty. In all other 
cases, infinitely big and rather unwieldy matrices have to be handled in order to find the results 
that interest us. These matrices furthermore depend on the one-point functions, which have to be 
solved for — a task that has been successfully completed only in a very limited number of settings 
[5] , none of them resembling the real world. 

On the other hand, the formalism has some as yet unexplored strength. Before writing down 
any equation, one must choose at what sites to take the functional derivative. This has no 
equivalent in ordinary Schwinger-Dyson equations and is due to the field variables in the path 
integral differing from the one we compute correlation functions of. Due to this choice, one has the 
liberty to search for different recastings of the equations, and it may very well be that a solution 
would come within reach if only the right choice is made. This will be left for further research. 

Another important avenue of further research would be the inclusion of quarks. Due to their 
fermionic nature, however, it is less natural to treat quarks on the lattice — due to the fermion 
doubling problem. Several techniques exist to avoid this problem, but always at the cost of 
increased complexity. 

Acknowledgements 

This work is supported by the Generalitat Valenciana under grant Prometeo/2008/004 and by the 
Spanish MICINN under grant FPA20 11-23596. 

References 

[1] Y. Nambu. QCD and the string model. Phys.Lett, B80:372, 1979. 

[2] E. Corrigan and B. Hasslacher. A functional equation for exponential loop integrals in gauge 
theories. Phys.Lett, B81:181, 1979. 

[3] J.-L. Gervais and A. Neveu. The quantum dual string wave functional in Yang-Mills theories. 
Phys.Lett., B80:255, 1979. 

[4] Yu.M. Makeenko and A. A. Migdal. Quantum chromodynamics as dynamics of loops. 
Nucl.Phys., B188:269, 1981. 

[5] T. Eguchi. Strings in \J{N) lattice gauge theory. Phys.Lett, B87:91, 1979. 

[6] D. Weingarten. String equations for lattice gauge theories with quarks. Phys.Lett., B87:97, 
1979. 

[7] D. Focrster. Yang-Mills theory — a string theory in disguise. Phys.Lett., B87:87, 1979. 

[8] D. Friedan. Some nonabelian toy models in the large N limit . Commun.Math.Phys., 78:353, 
1981. 

[9] A. A. Migdal. Loop equations and 1/N expansion. Phys.Rept, 102:199, 1983. 



15 



[10] W. Kerler. Schwinger-Dyson equations and currents in lattice gauge theory. Phys.Lett., 
B100:267, 1981. 

[11] D.J. Gross and E. Witten. Possible third-order phase transition in the large- iV lattice gauge 
theory. Phys.Rev., D21:446, 1980. 

[12] G. PafFuti and P. Rossi. A solution of Wilson's loop equation in lattice QCD2. Phys.Lett., 
B92:321, 1980. 



16 



